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ABSTRACT 

We consider the effects induced by the presence of hot and cold spots on the source star in the 
light curves of simulated microlensing events due to either single or binary lenses taking into 
account the rotation of the source star and the orbital motion of the lens system. Our goal is 
to study the anomalies induced by these effects on simulated microlensing light curves. 
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1 INTRODUCTION 

Gravitational microlensing has proven to be an exceptional tool to 
investigate several astrophysical phenomena. It is most often used 
to probe the lens system (located between the source and the ob¬ 
server) that may be a star, possibly with its planetary system, a 
MACHO (Massive Astrophysical Compact Halo Object) or even a 
free-floating planet. Moreover, gravitational microlensing behaves 
as a powerful natural telescope that allows to investigate the source 
star located at kpc distance from Earth (see, e.g., Gould 2001). The 
source star can be best studied when it crosses the caustics and 
moreover large finite source effects are present, that is when the 
projection of the source disc on to the lens plane is sufficiently 
large. This happens when the lens and/or the source are close 
enough or the source is a red giant star. 

The main applications of gravitational microlensing to stellar 
physics are the study of the limb-darkening of source stars (Witt 
1995) and of irregularities, like cold and hot spots, on their sur¬ 
face. Heyrovsky & Sasselov (2000) and Hendry et al. (2002) stud¬ 
ied the lensing of a spotted star by a single lens, Han et al. (2000) 
and Chang & Han (2002) later extended the study to the binary lens 
case. However, none of these papers has taken into account the pos¬ 
sibility that the source star and/or the binary lens system rotate. The 
aim of the present work is that of studying the anomalies induced 
by these effects on simulated microlensing light curves. 


2 EVENT SIMULATION 

In exoplanet search it is well known that the presence of stellar 
spots can induce features either in the observed lightcurves and in 
the radial velocity profiles, and may mimic the signal due to one 
or more planets orbiting around the star (Queloz et al. 2001). In¬ 
deed, recently that effect generated a false identification of a planet 
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in the initial analysis of the microlensing event MOA-2010-BLG- 
523 (Gould et al. 2013), as due to the intrinsic variability of the 
source star. 

This has pushed us to consider in detail the effects on simu¬ 
lated microlensing events of one or more starspots on the surface 
of the source star. 

As usual in describing microlensing events, we use use the 
Einstein radius Re and the Einstein time fg as the typical scales for 
lengths and times. They are defined as follows: 


Re - 


4GM D h D, 


D s 



v± 


(1) 

( 2 ) 


where M is the total mass of the lens system, D\ is the distance from 
the observer to the lens, D s is the distance from the observer to the 
source star, D ls = D s - D\ is the distance between the lens system 
and the source star, v ± is the speed of the source relative to the lens, 
perpendicular to the line of sight. In the following, all lengths and 
times are expressed in units of Einstein radius and Einstein time 
respectively. 

We calculate the amplification A s of the spotted star as the 
weighted average of the amplification A(r ) over the star disc S , 
using as weight the surface brightness /(r), which will be defined 
below, 


/ 5 A(r)/(r)d 2 , 

jT/(r)d 2 


The amplification A(r) can be either the amplification for a single 
lens (see, e.g., Paczynski 1986), 


A(u) = 


u 2 + 2 
u Vm 2 + 4 


(4) 


being u the projected distance between the lens and the source in 
units of Einstein radius, or the amplification induced by a binary 
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lens (Mao 2012). In the latter case, we calculate the amplification 
of the source by combining two different methods: 

• inverse ray-shooting (Schneider & WeiB 1986; Kayser et al. 
1986) when the source is close to the caustics; 

• Witt & Mao (1995) method when the source is far enough 
from the caustics. In particular, in order to solve the fifth- 
order polynomial equation we employed the root finder 
by Skowron & Gould (2012). 


The surface brightness profile f(r) in equation (3) is defined as 
I l(r) outside the spot 

fir) = { F (5) 

I c l(r) inside the spot 

where l(r) is the brightness of the unspotted star, and c is the con¬ 
trast parameter, that is the ratio between the brightness of the spot 
and the unspotted surface. The case c > 1 corresponds to a hot spot, 
c < 1 is for cold spots. To model the source star surface brightness, 
we adopt the linear limb-darkening law (Afonso et al. 2000) 


Hr) = Hr) = 


1 -T 




( 6 ) 


where r is the distance from the centre of the star disc, T,i is the 
limb-darkening coefficient at the wavelength A, p is the radius of 
the source, and F = F/np 1 is the total flux F of the star averaged 
over the stellar surface. The star spot modellization is described in 
appendix A. Integration over the star disc is carried out by using the 
Cuba library, developed by Hahn (2005). Since the stellar spots can 
be quite small, we have to be sure that the region around each spot is 
well sampled. Thus we chose in particular the Divonne algorithm, 
since it can allow us to increase sampling around the position of 
given points of the integrand function, the starspots in our case. 

Differently from previous works, in this paper we consider 
also the orbital motion of the lens system and the rotation of the 
source star. In the next section, we present the results of the simu¬ 
lated microlensing events to highlight the effects of source spinning 
and binary lens orbital motion. 


3 SIMULATIONS 
3.1 Single lens 

In Fig. 1 we show the results of a simulation for the single lens 
case. We set the ratio D\/D s between the distance to the lens and to 
the source to be 0.5; the source projected radius p = 0.0075 and its 
rotation period is 0.61 t E . The distance of closest approach between 
the point-like lens and the centre of the source is u 0 = 0. The two 
hot spots on the source have radius f? spots = 0.03 p and contrast 
parameter c = 4. One of the spots has latitude = 0, and so the 
lens passes right over it; the other spot has latitude 8 2 = tt/ 21 . 

It should be noted that, in this case, even if the star had 
a shorter period, the probability to see other spot-induced peaks 
would be negligible. In particular, it would be zero if the lens does 
not move parallel to the equator of the source, since the projected 
distance between the lens and the spot must be as small as possi¬ 
ble for the spot-induced peak to be visible (see the lower panel of 
Fig. 1). 


3.2 Binary lens 

In Fig. 2 our simulation results for two binary lens cases are shown. 
For simplicity, we assumed the orbit of the lens systems to have ec¬ 




Days 


Figure 1 . Simulations of a microlensing event by a single lens. In the up¬ 
per panel the source, with its hot spots, and the trajectory of the lens are 
represented. In the middle panel it is shown the amplification of the source 
with and without the spots. In the lower panel the fractional deviation of the 
amplification from that of a spotless event (A s - Aq)/Aq and the projected 
distance between the lens and each spot are shown. 


centricity e = 0 and to be face-on. The distance of closest approach 
between the source and the centre of mass of the lenses is uo = 0 . 1 . 
The source trajectories are shown in the upper panels of Fig. 2. 

For comparison purposes, in the left-hand panel of Fig. 2 we 
choose the system parameters of the event simulated to be quite 
close to the binary events simulated by Han et al. (2000). The lens 
system is constituted by two equal mass objects, so their mass ratio 
is q = 1, and a projected separation of b = 1. The ratio D\/D s is 
equal to 0.5. The orbital period, when non null, is set to 10.1 t E . The 
source star has a projected radius of p = 0.1 and its rotation period 
is 1.23 ?h- The limb-darkening coefficient used in the simulations is 
Li = 0.45. The spot, with centre on the equator, is a cold spot with 
c = 0.1 and radius equal to 0.2 p. 

In the right-hand panel of Fig. 2, we consider a simulated 
event with more realistic parameters, in particular smaller spots 
and a larger contrast parameter c, similarly to the real stellar 
spots detected with the transit method (see, e.g., Mancini et al. 
2014; Wolter et al. 2009). In this case the residuals are necessar¬ 
ily smaller, but the effects are nevertheless interesting to be stud¬ 
ied. The lens system components have the same mass (q = 1), and 
projected separation b = 1. Here, the ratio D\/D s has been fixed 
to 7/8. When the orbital motion of the lenses is active, the orbital 
period is 6 t E . The source star has a rotation period of 0.15 t E , and a 
projected radius of p = 0.1. The limb-darkening coefficient used in 
the simulation is r,i = 0.45. In this simulations, there are five cold 
spots with c = 0.25 and radius equal to 0.05p. Three of the spots 
are on the same side of the source sphere, with latitude 6\ = 0 , 
#2 = -7t/12, d 3 = 7T/ 18, and initial longitude = 0, <p 2 = —7t/8, 
(p 3 = 7 t /12 respectively. The other two spots are on the opposite 
side, so they cannot be seen in Fig. 2, have latitude 64 = 65 = 0, 
and initial longitude </> 4 = 2n/3, and <p s = 47 t/ 3 . 
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Figure 2. Simulations of microlensing events by a binary system. We simulated events with and without orbital motion of the lenses and rotation of the source 
itself. In the upper panels, the amplification map and the trajectories of the source in the lens plane t; — r) are shown in the case of no orbital motion of the 
system (red straight line) and if the lens system rotates (gold line). In the middle panels there are the amplification curve for the different trajectories and in the 
lower panels the residuals between the case of a spotted source and the unspotted one are reported. In the left-hand panels, a source star with one coldspot on 
its surface is considered, while in the right-hand panels, we show the results for a source star with five coldspots (only three initially on the visible side). We 
refer to the text for the parameters used in these simulations. 


4 RESULTS AND DISCUSSION 

The crossing of a caustic by a spotted star in microlensing events 
introduces extra peaks or dips in the amplification curve, which 
may be confused with planetary signatures. The spot-induced fea¬ 
tures are always very close to the main peak, because the secondary 


peaks or dips must fall within the angular size of the source. As 
already noted by Han et al. (2000), the probability to detect stel¬ 
lar spots is higher in the binary lens case than for the single lens. 
So, we simulated microlensing events by single and binary lens 
systems of sources with hot and/or cold spots. We also took into 
account the spin of the source star and the orbital motion of the 


MNRAS 000, 1-5 (2015) 

































































4 M. Giordano et al. 


lens binary. Indeed, we emphasize that the probability of detecting 
starspots tends to increase if the source rotates substantially dur¬ 
ing the lensing events (especially in the binary lens case) since the 
caustics are extended. 

The repetition of the same feature in the blue and gold 
residuals of the simulation in the right-hand panel of Fig. 2 
could, in principle, make it possible to estimate the rotation 
period of the source with a generalized Lomb-Scargle peri- 
odogram (Zechmeister & Kiirster 2009).' For a binary lens system, 
this does not even require the source to move parallel to its equa¬ 
tor in the lens plane, due to the fact the caustic is extended. This 
technique is already used for transit lightcurves, as tracking the 
change in position of a starspot allows one to recover not only the 
orbital obliquity of the host star with respect to the orbital plane, 
but also to get strong constraints on the stellar rotation period (see, 
e.g., Mancini et al. 2014). 

A possible way to distinguish between a caustic crossing peak 
and a spot-induced peak (c > 1 ) or dip (c < 1 ) is to observe 
the microlensing event in two different bands. Indeed, the con¬ 
trast parameter c strongly depends on the frequency of observa¬ 
tion (Silva 2003) and the presence of a spot characterized by a 
contrast parameter value much different than unity noticeably af¬ 
fects the microlensing lightcurve, being this effect more evident at 
shorter wavelengths. Of course, we remark that observing such sig¬ 
natures requires high precision and high cadence photometry and, 
possibly, observing the same target with a network of dedicated 
telescopes (such as KMTNet, see, e.g., Park et al. 2012). 

Finally, we remind that the optical depth for binary lens fea¬ 
turing star spots is 


= /1/2/3T 


(7) 


where r 4.48 x 1CT 6 is the microlensing optical depth towards 
the Galactic bulge, /j ^ 0.184 is the fraction of giant stars in the 
Galactic bulge observed by OGLE-III (Wyrzykowski et al. 2015), 
/2 - 0.054 is the fraction of binary lenses in the OGLE-III data 
with caustic-crossing features (Wyrzykowski et al. 2015), and, as 
assumed by Sajadian (2015) (see also Berdyugina 2005), / 3 — 0.01 
is the fraction of giant stars with stellar spots induced by a magnetic 
field stronger than 100 G. Here we are assuming that the efficiency 
of detecting star spot features with residuals larger than 10~ 3 to be 
about unity, based on reachable photometric precision for ground- 
based microlensing observations. It turns out that r, ^ 5x 10 -10 . 
The number of events with star spots can be estimated to be 


, 7 t Tobs^bg 

S= 2 <r E > 


( 8 ) 


where Mb g - 400 • 10 6 is the number of monitored stars in the 
Galactic bulge by OGLE-IV (Udalski et al. 2015), T ohs = 365 d, 
(f E ) =! 24.6 d (Wyrzykowski et al. 2015). Therefore, the number of 
events with star spots features is about 4.7 per year. We remark that 
the extra peaks or dips are always within the typical crossing time 
of the caustic. 
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APPENDIX A: MODELLING STAR SPOTS 

In this section we briefly describe the adopted modellization of the 
star spots. We assume the spots to be perfectly circular on the sur¬ 
face of the star, but we also take into account the effect of apparent 
distortion when they move towards the star limb due to the star 
rotation. 

In the sky plane, let (r, if/) be the polar coordinates of a generic 
point on the star disc, with radius p, then the three-dimensional 
Cartesian coordinates of this point on the star surface are 

x = pyjl -dip) 2 , (Al) 

y = rcost/r, (A2) 

z = r sin if/, (A3) 

where the y — z plane coincides with the sky plane, and the jc-axis is 
directed towards the observer. One can also calculate the colatitude 
9 and longitude tp of the point on the star surface with the usual 
relations 


9 = arccos(z/p), (A4) 

ip = arctan (y/x). (A5) 

Let (6i, ipi) be the angular spherical coordinates of the centre of the 
fth spot, with radius p,-, on the star surface. We need a criterion to 
determine whether a point with coordinates (9, tp) on the star surface 
is inside a star spot. For each spot, the quantities 

9' = 9 - (6>, - n/2), (A 6 ) 

V\ = V~ Vi ( A7 ) 

are the angular spherical coordinates of the point ( 9 , tp) centred on 
the centre of the spot. If cos tp'. > 0 then the spot and the point are 
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on the same side. The projected two-dimensional Cartesian coordi¬ 
nates of the point in this frame of reference are 

y'i = p sin ff i sin </>', (A8) 

z'=p cost?'. (A9) 

Thus, if the condition 

^iy'f + (z') 2 <Pi (Aio) 

is satisfied the point is inside the spot. In this way, the effect of dis¬ 
tortion of the spots follows naturally. We implemented this control 
in our integration routine in order to set the contrast parameter. 

This paper has been typeset from a T[:X/T>TpX file prepared by the author. 
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